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ABSTRACT 

We present the X-ray source catalog in the Subam/XMM-Newton deep survey. A continuous area of 1.14 
deg 2 centered at R.A. = 02hl8m and Dec. = -05d is mapped by seven pointings with XMM-Newton covering 
the 0.2-10 keV band. From the combined images of the EPIC pn and MOS cameras, we detect 866, 1114, 
645, and 136 sources with sensitivity limits of 6 x 10" 16 , 8 x 10" 16 , 3 x 10" 15 , and 5 x 10" 15 erg cm" 2 s" 1 in 
the 0.5-2, 0.5^1.5, 2-10, and 4.5-10 keV bands, respectively, with detection likelihood >7 (corresponding 
to a confidence level of 99.91%). The catalog consists of 1245 sources in total including 32 extended-source 
candidates. The averaged log N - log S relations are in good agreement with previous results, bridging the flux 
range between Chandra deep surveys and brighter surveys. The log N - log S relations show significant spatial 
variation among pointings on a scale of 0.2 deg 2 . Analyzing the auto correlation function, we detect significant 
clustering signals from the 0.5-2 keV band sample, which can be fit with a power law form (6>/# c )~ a8 with 
a correlation length of 6 C = 5.9"^ g,, when the integral constraint term is included. In the 2-10 keV band, 
however, the clustering is not significant with a 90% upper limit of 6 C < 1 .5". 

Subject headings: catalogs — diffuse radiation — galaxies: active — X-rays: galaxies — X-rays: general 



1. INTRODUCTION 



The Subaiu/XMM -Newton Deep Sur vey (SXDS; 



Sekigu chi et al.l2008l) is. along with COSMOS (IScoville et alJ 
2007t iHasinger et all |2007), one of the largest multi- 



wavelength survey projects with an unprecedented combi- 
nation of depth and sky area over a contiguous region of 
> 1 deg 2 . The main aims of the SXDS are to make an 
accurate measurement of the global properties of the universe 
without being affected by cosmic variance and to reveal the 
evolution of the large scale structure. The SXDS consists 
of a wealth of multi-wavelength data taken by the most 
modern observing facilities; X-ray imaging/spectroscopic 
data in the 0.2-10 keV ba nd taken by t he Euro pean Photon 
Imaging Camera (EPIC; Striider et al. 2001; Turner et al. 
1200 ll) onboard XMM-Newton dJansen et all 1200 ll) . multi- 
color (B, V, R, i', z') dee p optical images by S uprime-Cam 
on the Subaru telescope dFurusawa et al.l 12008). deep near 
infrared maps (J, H, K) observed as the Ultra Deep Survey 
(UDS) in the United Kingdom Infrared Deep Sky Survey 
(UKIDSS; lLawrence et al.ll2007l) . mid- and far-infrared data 

1 Department of Astronomy, Kyoto University, Kyoto 606-8502, Japan 
- Department of Physics & Astronomy, University of Leicester, Leicester, 
LEI 7RH, UK 

3 Jodrell Bank Centre for Astrophysics, University of Manchester, Oxford 
Road, Manchester, M13 9PL, UK 

4 Subaru Telescope, National Astronomical Observatory of Japan, Hilo, 
HI, 96720, USA 

5 Astronomical Institute, Tohoku University, Sendai 980-8578, Japan 

6 Astrophysikalisches Institut Potsdam, An der Sternwarte 16, 14482 Pots- 
dam, Germany 

7 Instituto de Fisica de Cantabria (CSIC-UC) Avda. de los Castros, 39005 
Santander, Spain 

8 National Astronomical Observatory of Japan, 2-21-1 Osawa, Mitaka, 
Tokyo 181-8588, Japan 

9 Astrophysics Research Institute, Liverpool John Moores University, 
Twelve Quays House, Egerton Wharf, Birkenhead CH41 1LD, UK 

10 Max-Planck-Institut fur extraterrestrische Physik, Giessenbachstrasse 
85748 Garching, Germany 



(3.6-160 /im) taken with the Spitzer Space Observatory 11 , 
the submillimeter (850 ^m) map b y the SCUBA Half Degree 
Extragalactic Survey (SHADES: fMortier et al.ll2005l) 12 . and 
the deep radio imag e (1.4 GHz) by the Very Large Array 
(ISimpson et al.ll2006l) . The SXDS field is centered at R.A. 
= 02hl8m and Dec. = -05d and the total field of the Subaru 
Suprime-Cam and XMM-Newton EPIC images covers an area 
of 1.3 and 1.14 deg 2 , respectively. The overa ll survey design 
and de tails of each survey are summarized in Sek iguchi et alJ 
d2008h . 

The data of XMM-Newton constitute a major component 
of the SXDS project. X-ray surveys are a powerful tool 
to trace the cosmological evolution of active phenomena in 
the universe, including Active Galactic Nuclei (AGNs) and 
clusters/groups of galaxies. The main constituents of X- 
ray sources that make up the X-Ray Background (XRB) 
are A GNs (for a recent review, see Brandt & Hasinger 
2005). Their domin ant populations are "obscured" AGNs 
(ISetti & Woltjerlfl989l) . where the central engine is blocked 
by dust and/or gas in the line of sight. In particular, sur- 
veys by hard X-rays above 2 keV are the most efficient ap- 
proach to detect these obscured AGN populations of various 
luminosity classes with least bias, thanks to their strong pen- 
etrating power against the photo-electric absorption of mat- 
ter and small contamination from stars. In fact, the surface 
density of AGNs detec ted in the most sensitive X-ray surveys 
dAlexander et al. 2003) far exceeds that detected in the optical 
bands (e.g.. |Wolf et al.l2003l) . thus providing us with the most 
complete and clean samples for AGN studies including heav- 
ily obscured, low luminosity ones that dominate the whole 
AGN populations by number. 

The tight correlation between the mass of a Super Mas- 
sive Black Hole (SMBH) in the galactic center and that of 
the galactic bulge dFerrarese & Merittl l2000t iGebhardt et al.l 

11 the Spitzer Legacy Survey and part of the Spitzer Wide-area InfraRed 
Extragalactic (SWIRE) survey 

12 The SXDS field is also the field of SCUBA-2 Cosmology Legacy Sur- 
vey. 
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2000) indicates strong links between the growth of SMBHs 
and star formation in the past. Since AGNs are phenomena 
that mark the process of growth of SMBHs in galactic cen- 
ters, to elucidate the evolution of AGNs is a fundamental issue 
for understanding the history of the universe. Combination of 
ultra-deep pencil beam surveys and large area surveys have 
revealed the cosmol ogical evolution of the X-ray luminosity 
function of AGNs dUeda et a ll 120031; lLa Franca et alj 120051: 
Bargeretal. 2005; Hasinger et al. 2005). The global accre- 
tion histor y of the universe is simi l ar to that of star form ation 
(see e.g., iFranceschini et al.1l 1999t. [Marconi et al. 2004]). An 
important finding is that the number density of more lumi- 
nous AGNs have a peak at higher redshifts compared with 
less luminous ones. This behavior is called "down-sizing" or 
"anti-hierarchical" evolution, which is opposite to a naive ex- 
pectation from the standard structure formation theory of the 
universe. Similar tr ends were also obtained in the star for - 
mation history (e.g. JCowie et al.lll996t iKodama et al.ll2004l) . 
These facts imply two modes of "co-evolution" of galaxies 
and SMBHs for different masses. It could be explained by 
the fe edback from supernova and AGNs (e.g., iGranato et alj 
2004). Deep multi-wavelength surveys enable us to observe 
various populations in different evolutionary stage, giving the 
best opportunities to investigate how galaxies and SMBHs co- 
evolved with cosmic time. 

Despite this major progress on AGN evolution in recent 
years, there are several critical issues to be resolved by current 
and future X-ray surveys, even if we limit the scope only to X- 
ray detected AGNs. Immediate objectives include (1) whether 
or how the fraction of obscured AGNs evolves with redshift, 
(2) the number density of Compton thick AGNs (those whose 
line-of-sight absorption exceeds A^h > 10 24 cm" 2 ), (3) the 
evolution of the AGN luminosity function at high redshifts 
(z 4) to be compared with the results from optical QSO sur- 
veys, and (4) the evolution of clustering properties of AGNs 
as a function of luminosity and type. The SXDS X-ray sur- 
vey, one of the few "wide" and "deep" surveys, will give im- 
portant steps to answer these questions and to establish the 
average properties of the AGN populations. Population syn- 
thesis models of the XRB utilize the results from deep surveys 
in a limited area at the faintest flux levels dUeda et al.ll2.003t 
iGiili et al] |2007). which would be subject to cosmic variance 
when we discuss the XRB spectrum with wl0% accuracy. 
This is crucial when we discuss a relatively minor contribu- 
tion of some X-ray sources, such as Compton thick AGNs. To 
constrain the statistical properties of rare objects, such as high 
redshift QSOs, much larger cosmic volume is necessary than 
in the currently available deep surveys. By surveying a contin- 
uous area with a large depth, we can also study the evolution 
of the large scale structure probed by AGNs. It helps us to 
understand the physical conditions that triggers AGN activity 
for different luminosity class, type (1 or 2), and redshift. At 
the flux limits of the SXDS, we can detect the most dominant 
populations of X-ray sources that contribute to the XRB. In 
all of these studies, high completeness of identification (red- 
shift determination) is crucial. The multi-wavelength data are 
particularly useful in identifying X-ray sources even if optical 
spectro scopy is difficult for a population of "optically faint" 
AGNs (lAlexanderet al.ll200l . 

In this paper, we describe the X-ray data of the SXDS based 
on the XMM-Newton survey performed in 2000 and 2003. We 
present the whole X-ray catalog, and basic statistical prop- 
erties of the detected X-ray sources. The X-ray data have 
been used in a number of studies in the SXDS projects, de- 



tection of ultra lum inous X-ray sources in nearby galaxies 
dWatson et alj |2005|). mid-IR and radio sele cted Compton- 
thick AGNs dMartfnez-S ansigre et al. 2007 ), radio sources 
dSimpson et al.l 120061) . Lyman a sources JSaito et al.l 120081; 
lOuchi et al. 2008 |), and optical variability selected AGNs 
(Morokum a et al.ll2008l) . A series of papers on optical iden- 
tification of the X-ray sources are forthcoming (Akiyama et 
al., in preparation). The paper is organized as follows. § 2 
summarizes the observation of the XMM-Newton. In § 3 we 
describe the detailed procedure of data analysis, and present 
the source list and its statistics. § 4 and § 5 shows the results 
of log N - log S relations and auto correlation functions, re- 
spectively. In § 6 long-term time variability of X-ray sources 
is studied in a selected field. § 7 summarizes the conclusions. 
Throughout the paper, the cosmological parameters (Ho, Sl m , 
n x ) = (100/z km s" 1 Mpc" 1 , 0.3, 0.7) are adopted. 

2. OBSERVATIONS 

The SXDS field, centered at RA. = 02hl8m and Dec. = - 
05d, was mapped with XMM-Newton in the 0.2-10 keV band 
by seven pointings, one deep (nominal exposure 100 ksec) ob- 
servation in the center surrounded by six shallower (50 ksec 
each) ones. XMM-Newton carries three EPICs, pn, MOS1, 
and MOS2, each has a field of view (FOV) of about 30' x 30'. 
The pointing positions are arranged so that the combined X- 
ray map, having a "flower petal" pattern, covers essentially 
the whole region of the optical images taken with the Sub- 
aru telescope, mosaic of 5 Suprime-Cam images in each band 
dFurusawa et alj 120081) . There are small overlapped regions 
between these pointings, making the X-ray image continu- 
ous without gaps over an area of about 1.14 deg 2 , although 
the sensitivity is not completely uniform over the entire map 
due to different exposures and instrumental effects, such as 
vignetting (see § 4.1 for details). 

Table 1 gives the observation log, including the pointing 
position, observation time, and net exposure obtained after 
screening out high background time. We designate the point- 
ings as SDS-1 for the central field and SDS-2 through SDS-7 
for the surrounding ones in a clockwise direction. The ob- 
servations were performed in 3 discrete epochs, the first from 
2000 July 31 to August 8 (for SDS-1, 2, 3, 4), the second from 
2002 August 8 to 12 (for SDS-5, 6, 7), and the third on 2003 
January 7 to supplement the unfulfilled exposure of SDS-4. 
The observation of each pointing were performed continu- 
ously except for short intervals (several hours) in SDS-1 and 
SDS-3 and a long interval (two and half years) in SDS-4. The 
"thin" filters were used for the three cameras in all the obser- 
vations. 

3. ANALYSIS 
3.1. Image Production 

The reduction of the EPIC data was performed by using 
the Science Analysis System (SAS) software. We used the 
SAS version 6.1.0 13 for image creation and all the subse- 
quent analysis, while an earlier version of SAS (that very 
similar to version 5.3.3) was utilized to produce the event 
files from which we start the image analysis, processed in the 
XMM-Newton Survey Science Center (SSC) Pipeline Process- 
ing System (PPS) facility in University of Leicester. 

As the first step, we created sky images with a bin size of 
4 arcsec in four energy bands, 0.3-0.5, 0.5-2, 2-4.5, 4.5-10 

13 The latest Current Calibration Files (CCF) as of 2005 January 31 were 
used. 
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TABLE 1 

LOG OF THE XMM-Newton OBSERVATIONS IN THE SXDS FIELD 



Field 


Obs. ID 


RA a 


Dec a 


PA b 


Start 


End 


Exposure 






(deg) 


(deg) 


(deg) 


(UT) 


(UT) 


(sec) 


SDS-1 


0112370101 


34.47838 


-4.98117 


73.6 


2000/07/31 22:38 


08/01 14:04 


40361 




0112371001 


34.47883 


-4.98050 


74.5 


2000/08/02 21:13 


08/03 13:15 


42473 


SDS-2 


0112370301 


34.87900 


-4.98047 


74.8 


2000/08/04 21:05 


08/05 14:00 


40220 


SDS-3 


0112370401 


34.67971 


-4.63422 


75.7 


2000/08/06 07:53 


08/06/17:38 


14341 




0112371501 


34.67958 


-4.63414 


75.7 


2000/08/06 20:58 


08/06/23:20 


4080 


SDS-4 


0112371701 


34.27946 


-4.63394 


76.0 


2000/08/08 07:57 


08/08 15:02 


21277 




0112372001 


34.32404 


-4.66858 


245.7 


2003/01/07 04:41 


01/07 11:56 


25940 


SDS-5 


0112370601 


34.07854 


-4.98164 


76.1 


2002/08/12 06:06 


08/12 18:52 


34377 


SDS-6 


0112370701 


34.27754 


-5.32797 


74.9 


2002/08/08 15:28 


08/09 04:43 


46802 


SDS-7 


0112370801 


34.67792 


-5.32856 


74.9 


2002/08/09 05:52 


08/09 19:02 


36803 



' d Mean pointing position (J2000) of the optical axis. 
b Mean position angle 

c Net exposure for the pn data after screening 

keV, separately for each observation ID and detector. The 7.9- 
8.3 keV band was ex cluded to avoid th e instrumental back- 
ground of Cu K-lines (Striider et al. 2001). With the standard 
event selection, we accumulated good photon events with 
"patterns" of through 12 (MOS1 and MOS2) or through 
4 (pn). The time region of background flare was excluded by 
examining the light curve. This leaves a net exposure of about 
60-80% of the allocated observing time, as listed in Table 1 . 

To achieve the best positional accuracy, we corrected the 
absolute astrometry of the EPIC image of each obs. ID with 
respect to the Subaru R-band image in the following way. We 
performed source detection from each observation, and cross 
correlated t he obtained (tentative ) source list with the opti- 
cal catalog ( Fur usawa et al.ll2008l) to calculate the mean posi- 
tional offsets. Then we corrected the original images for these 
offsets, which were found to be 0"- 1.3", consistent with the 
estimated error in the absolute astrometry of XMM-Newton. 
This correction also minimizes systematic errors caused by 
summing up multiple observations for the case of SDS-1, 3 
and 4. 

3.2. Source Detection 

We produce source lists separately from different pointings 
(SDS-1 through 7). Later, these are merged into one list by 
excluding overlapping regions to define a statistically inde- 
pendent sample from the whole data (§ 3.3). The overall 
flow of the source detection is similar to the one adopted in 
the SSC PPS used to produce the 2XMM catalog (Watson et 
al., in preparation), although several improvements are em- 
ployed here. The main differences are as follows; (1) we per- 
form source detection to the summed image of pn, MOS1, 
and MOS2, not to individual ones. (2) We optimize several 
parameters of the spline interpolation in making background 
maps. (3) Special care is paid to separate neighboring sources. 

3.2.1. Background Maps 

For each observation, energy band, and detector, we cre- 
ate a background map from the image by spline interpolation 
with the SAS program esplinemap after excluding sources 
found with eboxdetect, which perform simple source detec- 
tion based on a cell detection algorithm. After optimization of 
parameters, we verify that the obtained models well represent 
the profiles of the background (i.e., the unresolved XRB plus 
non X-ray background) over the whole FOV without large de- 
viation from the data. 

3.2.2. Summation of Images 



We sum up the pn, MOS1, MOS2 images and background 
maps for each energy band and each "pointing", which means 
also combining those of multiple observations for SDS-1, 3 
and 4. This enables us to achieve the best sensitivity from 
the whole available data in the simplest manner. At the same 
time, we can avoid technical problems in the position and flux 
determination that would be caused by image gaps between 
CCD chips in a single detector. 

To obtain exposure maps, we first calculated them for each 
instrument in narrow energy bands of every 0.5 keV, using the 
SAS task eexpmap. Then, we took their count-rate weighted 
average, normalized to the pn count rate in a given survey 
band. The count rate distribution is calculated through the 
energy response of each instrument assuming a power law 
photon index of 1.5 above 2 keV and 2.0 below 2 keV. The 
dependence on assumed photon indices is negligibly small in 
the analysis. The cross-calibration of effective area between 
the pn and MOS cameras is known to be accurate within 5% 
level, and its uncertainty does not affect the results. 

With this procedure, we are able to treat the summed im- 
ages and exposure maps as if they were created from a single 
pn detector regardless of the detected position. In the follow- 
ing analysis, we present the flux of a source in terms of a 
"pn-equivalent" count rate. Strictly speaking, this treatment 
may not be perfectly accurate in the image fitting process de- 
scribed below, because the positions of the optical axis in the 
combined image are not common for different detectors and 
observations, which could affect a precise modeling of the 
combined Point Spread Function (PSF) of the mirrors. To es- 
timate maximum systematic errors caused by this effect, we 
make the same analysis by changing the position of the op- 
tical axis to that of a MOS camera. We find that the fluxes 
of sources obtained in the two analyses match each other by 
1.3% level, confirming that this approximation is justified at 
this accuracy. 

Figure 1 shows the pseudo three-colored X-ray image of 
the SXDS field combined from the 7 pointings. This image 
is made by the export command in the IRAF package 14 . As 
the inputs, we use the exposure-corrected, background sub- 
tracted pn+MOSl+MOS2 images in the 0.5-2, 2^1.5, and 
4.5-10 keV bands, corresponding to red, green, and blue col- 
ors, respectively. They are smoothed by a Gaussian profile 
with a lcr width of 4 arcsec over the entire field. Larger color 

14 IRAF is distributed by the National Optical Astronomy Observatories, 
which is operated by the Association of Universities for Research in Astron- 
omy, Inc. under cooperative agreement with the National Science Foundation 
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fluctuation is apparently seen in outer regions compared with 
the central region due to the shorter exposure. 

3.2.3. Maximum Likelihood Fit 

The source list from each pointing is obtained by the maxi- 
mum likelihood fit applied to the detector co-added images in 
the four energy bands by using the SAS task emldetect. The 
emldetect program reads images, background maps, and ex- 
posure maps in a single or multiple energy bands, and makes 
a simultaneous fit to the images with a model consisting of 
PSFs over the background map, based on the input list of 
source candidates. It returns the fitted position, vignetting- 
corrected count rate, and detection likelihood (hereafter called 
"ML" standing for Maximum Likelihood) of each source. The 
obtained count rate corresponds to the total flux in the entire 
PSF. 

To make the list complete, this process is iterated twice 
in the following manner. We first perform cell detection by 
eboxdetect with a cell size of 20" x 20" to obtain a list of 
source candidates with minimum likelihood of 3.2, which is 
used as an input to the first run of emldetect. In both pro- 
grams we allow the flux and position of a source to be free 
parameters. We find that in some cases the above cell size 
is too large to completely detect sources in a high density re- 
gion. To supplement this, we also perform cell-detection with 
a 12" x 12"cell size, and add new source candidates detected 
with high significance into the source list obtained above. In 
this stage, by human inspection, we pay attention not to in- 
clude obviously fake sources such as a part of a PSF tail from 
nearby bright sources, diffuse emission, and false detection 
close to gaps between CCD chips. The merged list is again 
input to emldetect, producing a final source list. 

Through the fitting process, we find that there are sources 
that are likely to be confused by the PSF of neighboring ones 
in crowded regions. To measure the fluxes of these sources 
with the best accuracy, we perform an image fit in each region 
by allowing the positions and fluxes of the multiple sources 
(up to three) to vary simultaneously in order to solve the cou- 
pling of the fluxes self-consistently. Their angular separation 
is typically 12"- 24". Their fluxes are replaced with the new 
values obtained here, and the ML values are calculated for 
each source at the fixed position. The number of pairs of such 
possibly confused sources are about 10 per pointing, and are 
marked by flag "C" in the last column of the source catalog 
(Table 2). 

The statistic that emldetec t uses to fit th e source parameters 
is the C statistic defined by ICash I d!979f) . For a Poissonian 
probability distribution, appropriate to the measurement of X- 
ray events by XMM-Newton, Cash's statistic takes the form 

N M 

C = 2 y~^ (gi, j — njj In Cjj) + a constant, 

i=l 7=1 

where ey is the expected number of X-ray events in the ith 
pixel and jth energy band, n,-j is the detected number of 
events in that pixel/band, and the sum is over the N pixels 
within the detection region and the M energy bands used (here 
= 4 as stated). The value of e/ j is obtained by adding the value 
of the source model at that pixel/band to the expected con- 
tribution from the background; ey is thus a function of the 
source model parameters. 

After arriving at those values of the source parameters 
which minimize C, the detection likelihood (formally, the 
probability of the null hypothesis) for those optimum values 



is then calculated. Cash's prescription for this is to form the 
difference 

AC = C nu ll — Cbest 

where C nu u is C calculated by setting the amplitude of the 
source model to zero and Cb es t is the minimum result returned 
by the fitting routine. According to Cash's theory, AC is 
distributed approximately as \ 2 f° r v degrees of freedom, 
where v is the number of fitted parameters. The probabil- 
ity p(x 2 > AC) of obtaining the calculated value of AC or 
greater by chance fluctuations of the detected background can 
therefore be obtained in terms of the complementary incom- 
plete gamma function Q as follows: 

p( X 2 > AC) = 2W2,AC/2). 

The value of ML is set within emldetect to equal —ln(p). Note 
that the statistics must be treated with caution when the total 
number of photon counts used for the fit are very small ( 
9), which is not the case in our analysis. 

The ML value as calculated in the above way is viewed as 
more sensitive than the simple box-detection statistic used by 
eboxdetect, because it uses information about the source PSF 
to help reject random fluctuations of the detected background. 
Offsetting this is a possible bias introduced by the fact that 
emldetect is not given 'free rein' over the whole field, but is 
applied only to a small number of restricted areas of the field, 
at locations of candidate sources already found by eboxde- 
tect. If the detection threshold of eboxdetect is set to too high 
a value, there is a risk that a real source with a good shape, 
but not many counts, will be missed by the combined detec- 
tion process. In order to avoid this bias, one ought therefore 
to run the preceding eboxdetect with a detection threshold set 
deliberately low. On the other hand, since emldetect must 
go though the computationally-intensive fitting process with 
each of its candidates, it will be impractical to provide it with 
too long a list of candidates. The respective values of 3.2 and 
7 adopted here for the eboxdetect and emldetect likelihood 
cutoffs were chosen with this necessary compromise in mind. 

An additional advantage of using emldetect to determine 
the final source parameters is its ability to add the models of 
already fitted sources to an internal background map. Since 
the sources are processed in the order of their brightness, it 
is possible to take into account the background introduced by 
bright sources when fitting the fainter sources. 

3.2.4. Examining Source Extent 

Up to this stage, all sources are assumed to be point like. To 
constrain the spatial extent of the detected sources, we per- 
form an image fit by allowing the source extent to be a free 
parameter assuming a Gaussian profile. The source positions 
are fixed at the input ones. This procedure yields a list of ex- 
tended source candidates from the SXDS field (Table 3). For 
these sources, the fluxes listed in Table 2 correspond to those 
obtained in this process with consideration of the source ex- 
tent. However, as our procedure is essentially dedicated to de- 
tection of point like sources in the earlier stage, this list must 
be regarded to be incomplete. A more extensive approach 
to search for extended sources will be presented elsewhere 
(Finoguenov et al., in preparation). 

3.3. Source List 

Thus, we obtain seven source lists separately from differ- 
ent pointings (SDS-1 through 7), containing the information 
of position and fluxes with the ML values in the four energy 
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FIG. 1 . — The smoothed, 3-colored X-ray image of the whole SXDS field obtained from the energy bands of 0.5-2 keV (red), 2^.5 keV (green), and 4.5-10 
keV(blue). 



bands, 0.3-0.5, 0.5-2, 2-4.5, and 4.5-10 keV. In these tenta- 
tive lists, we keep all sources whose summed ML from the 
four bands is equal to or exceeds 5, which will be further 
screened in the following way. 

To make scientific research in comparison with previous 
studies, it is quite useful to define a sample selected in stan- 
dard bands such as 0.5-4.5 keV (hereafter "XID band") and 
2-10 keV band ("hard band"). Although the combined flux 
and ML can be calculated from those in the two individual 
bands (i.e., 0.5-2 and 2-4.5 keV for the XID band, and 2- 
4.5 and 4.5-10 keV for the hard band), its statistical treat- 
ment would become quite complex; for instance, the sensitiv- 
ity limit cannot be uniquely determined as a function of ML at 
a given position, depending on the source spectrum. Hence, to 
obtain a well defined sample selected in the XID or hard band 
in the same way as in the four narrow bands, we calculate their 
ML values by repeating the likelihood fit using the 0.5^.5 or 
2-10 keV image with the background map summed from the 
two narrow bands. In the fit, the positions of the sources are 



fixed 15 and no new sources are considered. 

In this paper we adopt a threshold for the detection likeli- 
hood of 7 in a single band. Thus, from the tentative source 
lists, we select only those detected with ML>7 either in 
the 0.3-0.5 keV (ultra-soft), 0.5-2 keV (soft), 2-4.5 keV 
(medium), 4.5-10 keV (ultra-hard), 0.5-4.5 keV (XID), or 
2-10 keV (hard) band, to be included in the final list. Fur- 
ther, for statistical analysis using this sample, it is convenient 
to have a single source list merged from the seven pointings; 
as mentioned above, there are overlapping regions between 
different pointings in the outer edge of the FOV, where same 
sources are repeatedly detected. To exclude such duplication, 
we only refer to the results of a single pointing that achieves 
the highest sensitivity at a given position, based on sensitivity 
maps in the XID band (see below). 

15 Since we fix the position in the XID and hard bands, the corresponding 
AC value is smaller than those obtained in the four individual bands for a 
given ML, because the degree of freedom v is reduced by 2. This is also 
the case for the "c onfuse d" sources whose ML values are obtained by fixing 
the position (see jj |3.2.3t . We properly take these effects into account in the 
analysis of log N - log S relations and auto correlation function. 
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TABLE 3 
The list of Candidates of 
Extended Sources in the SXDS 



No. 


Extent Likelihood 


Extent 11 
(") 






Z.OZH U.J 


0051 


11.4 


17 9i 1 2 


0140 


6.8 


22. 0± 3.3 


0153 


9.2 


4 7^ o.4 


0156 


7.4 


4 2i 9 


0239 


4.4 


3 9i 1 


0280 


4.5 


19. 3± 3.3 


0285 


36.6 


260± 1.8 


0287 


4.9 


4.8± 1.3 


0292 


5.5 


4.1± 1.3 


0396 


4.0 


5.9± 1.8 


0441 


12.4 


12. 7± 1.4 


0453 


7.2 


12. 4± 2.2 


0514 


30.9 


2i q± 2 


0552 


6.5 


7_lj- 14 


0621 


6.5 


6.6± 1.4 


0622 


43 


4.6± 1.4 


0625 


5.7 


2.9± 1.8 


0646 


9. 1 


3 5i 8 


0647 


10.3 


9.4± 1.6 


0712 


6.6 


4.3± 0.8 


0784 


14.3 


33.6± 3.3 


0796 


20.5 


3.0± 0.3 


0829 


48.6 


11. 5± 0.6 


0852 


4.7 


3.3± 0.2 


0876 


19.6 


17.4± 1.6 


0889 


5.7 


17.2± 3.1 


0934 


10.0 


4.7± 0.4 


1152 


5.0 


U.8± 3.9 


1168 


4.6 


18.4± 2.8 


1169 


5.2 


3.7± 1.1 


1176 


31.4 


15.3± 1.1 



a The source extent when a Gaussian pro- 
file is adopted (lex). 



Table 2 gives the source list in the SXDS field complied 
in this way, sorted by right ascension (RA) and declina- 
tion (Dec): (col. 1) source number; (cols. 2 and 3) the X- 
ray source position as determined by XMM-Newton (RA and 
Dec); (col. 4) the statistical error in the position estimated 
through the maximum likelihood fit (root sum square of the 
lcr error in each direction); (cols. 5 through 10) the ML value 
in each energy band; (cols. 11 through 14) the vignetting- 
corrected count rate in each energy band; (cols. 15 trough 
17) the hardness ratios defined as HRl = (S-US)/(S+US), 
HRl = (M-S)/(M+S), HR3 = (U H - M)/(UH +M), where 
US, S, M, and UH are the count rates in the 0.3-0.5, 0.5-2, 
2-4.5, 4.5-10 keV band, respectively; (col. 18) the pointing 
ID at which the source is detected; (col. 19) the offset an- 
gle from the mean optical axis in the corresponding pointing; 
(col. 20) the total pn-equivalent exposure (sum of pn, MOS1, 
and MOS2) at the source position, corrected for the vignetting 
in the 0.5-4.5 keV band; (col. 21) the background rate in the 
0.5-4.5 keV band at the source position; (col. 22) the flag if the 
flux is determined by multiple-source fit to take into account 
possible source confusion with nearby sources (see § 3.2.3). 
The exposure and background rate (cols. 20 and 21) are aver- 
aged over the source extraction region with a fixed encircled 
energy fraction (68%) of the PSF. From the results of opti- 
cal identification (Akiyama et al., in preparation), we confirm 
that the positional errors given in this list are reasonable. In 
Table 3 we list the spatial extent and its likelihood of all the 
extended sources contained in Table 2. 

We detect 1245 sources in total. The numbers of detected 



sources in a single or two bands in any combination from the 
six energy bands are summarized in Table 4. According to 
simulation study 16 , the number of fake sources due to statis- 
tical fluctuation is estimated to be about 1.5 per pointing in 
each energy band at our likelihood threshold (ML>7), hence 

50 in the sum sample. Their fraction in the total number 
of detected sources depends on the selection band; it becomes 
the smallest (0.6%) in the 0.5-4.5 keV but largest (5%) in the 
4.5-10 keV band because of different sensitivities. 

We calculate the conversion factors from a count rate into 
a flux in each energy band by using the energy response 
of the pn. Galactic abso rption of Nu = 2.5 x 10 20 cm" 2 
dDickev & Lockmanll 19901) is taken into account. The values 
are summarized in Table 5 for a power law spectrum with var- 
ious photon indices (V). Figure 2(a) shows their dependence 
on the photon index assuming no absorption, whereas in Fig- 
ure 2(b) we change the absorption column density at zero red- 
shift for a fixed photon index of 1.8. In the same figures we 
also plot the hardness ratios HRl, HRl, and HR3 as a func- 
tion of spectral parameters. For the flux conversion, hereafter 
we assume T = 1.5 (with no absorption) for the 2^4.5 keV, 
4.5-10 keV, 0.5^1.5 keV, and 2-10 keV band, and T = 2.0 for 
the 0.3-0.5 and 0.5-2 keV band, unless otherwise mentioned. 

3.4. Spectral Properties of the Sources 

In this subsection, we summarize spectral properties of the 
SXDS sources based on hardness ratio analysis using the 
above source list. Figures 3(a) shows the flux distribution of 
sources detected in the 0.5^1.5 and 2-10 keV bands. As de- 
scribed below, our sample covers X-ray sources that constitute 
a major fraction of the XRB below 10 keV. Figure 3(b) shows 
the histogram of hardness ratio HRl and HR3 for those de- 
tected in the 0.5^4.5 and 2-10 keV bands, respectively. The 
flux versus hardness ratio plots are given in Figures 4 (0.5- 
4.5 keV flux versus HRl) and 5 (2-10 keV flux versus HR3). 
For clarity we do not attach error bars to each point in these 
figures; the typical lcr errors in HRl and HR3 are 0.22 and 
0.15 at flux levels of (2-4) x 10" 15 erg cm" 2 s" 1 (0.5^.5 
keV) and (1-2) x 10" 14 erg cm" 2 s" 1 (2-10 keV), respec- 
tively. To examine the mean spectral properties as a function 
of flux, we produce flux-weighted spectra in the four bands 
from flux-range limited samples, by summing up the count 
rate of the individual sources in each energy band. Table 6 
summarizes their HRl and HR3 values with corresponding 
power-law photon indices. We find evidence that the spec- 
tra of the X-ray sources becomes harder at fainter fluxes with 
a confidence level of 99.99% (HRl) and 96% (HR3) for the 
0.5-4.5 keV sample but with a lower confidence level of 84% 
(HRl) and 90% (HR3) for the 2-10 keV sample. Their mean 
slope is similar or even harder than that of the XRB, T ~ 1 .4 
over the 0.5-10 keV band at these flux levels. This result is 
consistent with the fact that we already resolved a major parts 
of the XRB and that at brighter fluxes where 20-30% of the 
XRB is produce d, the mean spec tra are much softer than the 
XRB spectrum dUeda et al .11 1999b . 

Figure 6 shows the color-color plot between HRl and HR3 
by using 208 sources detected in the 0.5^1.5 and/or 2-10 band 
whose statistical errors in both hardness ratios are smaller 
than 0.2. The solid curve tracks the color of a power law 
spectrum with varying photon indices from to 4, while the 
dashed curve does that with varying absorption column den- 
sities at z = 1 (a typical redshift of our sample; see Figure 14) 

16 see http://xmmssc-www.star.le.ac.Uk/newpages/lXMM/fig_5.3.html 
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TABLE 4 
Number of Detected Sources 



Detection Band 


0.3-0.5 0.5-2 


2^4.5 


4.5-10 


0.5^.5 


2-10 


One Band Only a 


0.3-0.5 keV 


357 339 


235 


73 


341 


253 


16 


0.5-2 keV 


866 


412 


113 


853 


466 


13 


2-4.5 keV 




487 


121 


474 


461 


4 


4.5-10 keV 






136 


125 


133 


3 


0.5-4.5 keV 








1114 


550 


169 


2-10 keV 










645 


78 



NOTE. — Number of sources commonly detected in the combination of two bands (or a 
single band) are listed. 

11 Number of sources detected only in one band given in the first column 



TABLE 5 
Flux Conversion Factors 



Assumed Spectrum 
Photon Index 


0.3-0.5 keV 


0.5-2 keV 


Flux Conversion Factors" 
2-4.5 keV 4.5-10 keV b 


0.5^1.5 keV 


2-10 keV b 


1.0 


0.164 


0.171 


0.574 


1.718 


0.305 


1.057 


1.1 


0.167 


0.171 


0.571 


1.696 


0.295 


1.026 


1.2 


0.171 


0.170 


0.568 


1.675 


0.285 


0.996 


1.3 


0.174 


0.170 


0.564 


1.653 


0.276 


0.967 


1.4 


0.177 


0.169 


0.561 


1.632 


0.267 


0.939 


1.5 


0.181 


0.168 


0.558 


1.612 


0.259 


0.912 


1.6 


0.184 


0.168 


0.555 


1.593 


0.250 


0.886 


1.7 


0.187 


0.167 


0.551 


1.573 


0.243 


0.862 


1.8 


0.190 


0.167 


0.548 


1.554 


0.236 


0.838 


1.9 


0.194 


0.167 


0.545 


1.536 


0.229 


0.816 


2.0 


0.197 


0.166 


0.541 


1.517 


0.223 


0.794 



a The conversion factor from the vignetting-corrected, pn-equivalent count rate into energy flux between the 
same band in units of [10~ 14 erg cirT 2 s~'] / [count ksec~']. The fluxes are corrected for Galactic absorption of 
N H =2.5 X 10 20 cm" 2 . 

b The 7.9-8.3 keV band is excluded in the count rate but included in the flux. 




Photon Index Column Density (10 22 cm 2 ) 



FIG. 2. — (Lower): The conversion factors from the vignetting-corrected count rate into flux in the same band, calculated for different spectra. The unit is 
[10~ 14 erg cirT 2 s _1 / (c/ksec)]. The solid, dot-dashed, long-dashed, short-dashed curves correspond to the 0.3-0.5 keV, 0.5-2 keV, 2-4.5 keV, and 4.5-10 keV 
bands, respectively. (Upper): The hardness ratio expected from the assumed spectra. The solid, long-dashed, and short-dashed curves correspond to HR\, HR2, 
and HR3, respectively. 

(a) Left: the parameter is a photon index of a power law with no absorption, (b) Right: the parameter is an absorption column density for a power law spectrum 
with a photon index of 1.8. 
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TABLE 6 

Average Hardness Ratio for Flux Limited Samples 





Sample 






Averange Spectra a 




Detection 


Flux Range 


Number of 


0.5^1.5 keV 


2-10 keV 


Band 


(Ur 14 ergs" 1 cm" 2 ) 


Sources 


<HR2> 


< To.5-4.5 > 


<HR3 > 


< r 2 -io > 


0.5^.5 keV 


0.5-2.0 


406 


-0.55 ±0.01 


1.55 ±0.02 


-0.27 ±0.02 


1.34 ±0.05 


0.5-4.5 keV 


2.0-8.0 


107 


-0.61 ±0.01 


1.71 ±0.02 


-0.33 ±0.02 


1.52±0.05 


2-10 keV 


1.0-4.0 


334 


-0.53 ±0.01 


1.50 ±0.02 


-0.23 ±0.02 


1.23 ±0.04 


2-10 keV 


4.0-16 


38 


-0.55 ±0.01 


1.55 ±0.02 


-0.29 ±0.03 


1.40 ±0.08 



a The averaged hardness ratio is derived the summed spectrum of the flux-limited sample. The corresponding power 
law photon index in the same band is also shown, which is corrected for Galactic absorption of Nu = 2.5 X 10 20 cm -2 . 
The error is la. 




FIG. 3. — Left: histogram of flux for the 0.5^1.5 keV (solid) and 2-10 keV (dashed) selected sample in the SXDS XMM-Newton catalog. Right: that of 
hardness ratio HR2 (solid) and HR3 (dashed) for the 0.5^1.5 keV and 2-10 keV selected sample, respectively. 
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-1.0 I 

0.05 0.10 



0.50 1.00 5.00 10.00 

0.5-4.5 keV Flux (10~ 14 ergs cm" 2 s" 1 ) 



FIG. 4. — The 0.5^.5 keV flux vs ffl?2 plot for the sources detected in 
the 0.5^.5 keV band. HR2 is defined as (M-S)/(M + S) where M and 5 is 
the vignetting-corrected pn-equivalent count rate in the 2^.5 keV and 0.5-2 
keV, respectively. The dashed line corresponds to a power law photon index 
of 1.8. 



-1.0 






0.5 1.0 2 5.0 10.0 20.0 

2-10 keV Flux (10~ 14 ergs cm- 2 s" 1 ) 



FIG. 5. — The 2-10 keV flux vs HR3 plot for the sources detected in the 
2-10 keV band. HR3 is defined as (H-M)/{H + M) where H and M is the 
vignetting-corrected pn-equivalent count rate in the 4.5-10 keV and 2^4.5 
keV, respectively. The dashed line corresponds to a power law photon index 
of 1.8. 

from Log A^h = to 24 for a fixed photon index of 1.8. As 
noticed, the spectra of some sources are not simply repre- 
sented by an absorbed power law model. In particular, there 
is a population of sources that are hard in the 2-10 keV band 
(e.g., HR3 > 0) but soft in the 0.5-4.5 keV band (HR2 < 0), 
indicating the presence of additional soft components. De- 
tailed spectral analysis using information of redshift will be 
presented in forthcoming papers. 

4. Log N - Log S Relation 

4.1. Sensitivity Map 

To make statistical analysis using the source list, such as 
determination of log N - log S relation, it is crucial to have 
reliable sensitivity maps for each pointing and energy band, 
i.e., we need to know a flux (count rate) limit as a function of 
position at a given detection criterion (i.e., ML > 7). Since 
we utilize the maximum likelihood algorithm, however, it is 
not trivial to calculate sensitivities by employing an analytic 
formula, unlike in the case of cell detection. Ideally, only a 
detailed simulation can give the correct estimate of a sensi- 
tivity at every position, since it depends on the background, 
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HR2 (2-4.5 keV vs 0.5-2 keV) 

FIG. 6. — The color-color plot between HR2 and HR3. Only sources with 
a lcr uncertainty less than 0.2 in both HR2 and HR3 are included in the plot. 
The solid curve (red) tracks the colors of an absorbed power law spectrum for 
a fixed photon index of 1 .8 with varying column densities (from log Nu cirT 2 
= 0, 21, 22, 23, 24 at z = 1), while the dashed one (blue) for varying photon 
indices from to 4 with no absorption. 



exposure map, and PSF in a complex way. To save comput- 
ing time, we here take an empirical approach by utilizing both 
simulation and analytical calculation, as described below. 

First, we directly estimate the flux limit by simulation at 
2-dimensional grid points with 40" spacing over the FOV. 
For a given position, we produce a simulated image where 
a point source is superimposed on an input background map, 
using the emldetect program. Since we are looking simply 
for a relation between the emldetect ML value and the flux 
of the sources, we do not taken into account Poisson count- 
ing noise in our simulation. We then perform a likelihood 
fit with emldetect to the simulated image and derive the ML 
value for the input source. Repeating the simulation assuming 
several different fluxes for the point source, we obtain an em- 
pirical relation between the input flux and output ML value. 
This relation enables us to calculate the sensitivity for a given 
threshold of ML at that position. 

Next, we produce full resolution (4" pixel) sensitivity maps 
by interpolating between the grid points in the following 
manner. The flux limit could change quite sensitively with 
the position, being affected by data gaps between CCD gaps 
or hot pixel regions. A simple interpolation of the flux limit 
based on the rough position sampling in the above simulation 
may not be sufficiently accurate. To estimate the precise po- 
sition dependence, we utilize analytical sensitivity maps cal- 
culated by the esensmap program, which is based on the cell 
detection algorithm. Generally, the detection likelihood ob- 
tained by emldetect (ML) differs from that defined in the cell 
detection (ML ce n) for the same flux. The ratio between ML 
and ML ce ii should be a function of position, primarily depend- 
ing on the size and shape of the PSF. Hence, we calculate the 
ML ce u value corresponding to the flux limit at each grid point, 
and interpolate them between the grid points as a linear func- 
tion of position. We finally convert ML ce n into a "flux limit" 
at each position, by referring to a set of sensitivity maps cal- 
culated by esensmap with different thresholds of ML ce n, thus 
producing the final sensitivity maps. To evaluate systematic 
errors in the sensitivity map, we compare the results when the 
positions of the grid points are shifted by 20" in each direc- 
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TABLE 7 
Best Fit Parameters for Log N ■ 



Logs Relations 



zzo.ooo ™d jsoto J-nao KkOT.O 

Right Aszeraion (JOOOj 

FIG. 7. — The sensitivity map in the 0.5—4. 5 keV band. The minimum flux 
threshold for the detection likelihood > 7 is given as a function of position in 
the sense that a darker color corresponds to a lower flux. SDS-2 is located to 
the east of the central field, SDS-1, and SDS-3 through 7 go in a clockwise 
direction after SDS-2. 

tion. We find that the difference in the flux limit is mostly 
within a few percent and <10% at maximum. We regard this 
(10%) as the maximum relative systematic error in the sen- 
sitivity map, which will be considered in the analysis of the 
spatial distribution of sources (§ I5.21 i. 

4.2. Results and Comparison with Other Work 

Superposing the sensitivity maps of the seven pointings in 
the common sky coordinates, we examine, position by posi- 
tion, in which pointing data the best sensitivity is achieved in 
the 0.5—4.5 keV band in overlapping regions. As mentioned 
above, we only refer to the results of a single pointing, at given 
sky coordinates, to obtain the combined source list and sensi- 
tivity maps to eliminate duplication. The resultant sensitivity 
map in the 0.5^4.5 keV band is shown in Figure 7. Figure 8 
shows the area curves il(S) for ML>7 (i.e., survey area where 
the detection completeness is guaranteed above the given flux) 
obtained from the combined SXDS field in the four energy 
bands, 0.5-2 keV, 0.5-4.5 keV, 2-10 keV, and 4.5-10 keV. 

Dividing the observed flux distribution by 0(5), we first ob- 
tain the log N - log S relations in the differential form, N(S). 
Here, we discard sources whose flux falls below the sensitiv- 
ity limit at the give position because we consider that such 
detection may not be reliable. Figure 9 shows N(S) in the 
four bands in units of number per square degree per 10~ 14 erg 
cm" 2 s" 1 , where the errors represent the la Poisson error in 
the number of sources in each flux bin. Following previous 
work, we fit them with a broken power law form, expressed 
as 



N(S) = 



f(|r r ' (s<s b) 



(i) 



The model has a break flux Sb, above and below which the 
slope is r„ and Tj, respectively. To estimate the parameters, 
we utilize a Maximum Likelihood method, so that the model 
best reproduces the observed flux distribution expected from 
the area curve. We fix T u at 2.5 in the 0.5-2 keV, 0.5-4.5 keV, 
and 2-10 keV bands, while in the 4.5-10 keV band, we adopt 



Band 


r u 


r d 


(i0- ]4 cgs) 


K 

(deg- 2 ) 


0.5-2 keV 
0.5-4.5 keV 

2-10 keV 
4.5-10 keV 


2.5 (fixed) 
2.5 (fixed) 
2.5 (fixed) 
z.oz_„ n 


1 63+0- 07 
1 63+ 006 
41^ 


i q 1+0.44 

LLy -0.1S 
1.0 (fixed) 


154+« 
177+52 

444+™ 
288+^0 



NOTE. — The error is lcr for a single parameter. 

a single power law by setting Yd = T u and 5b = 1 x 10~ 14 erg 
cm -2 s" 1 because the break is not evident in our flux range. 
The best-fit models are shown in solid lines in Figure 9. Ta- 
ble 7 summarizes the best-fit parameters, which are found to 

tent with pi 

ICarrera et al]|2007l:lBrunner et alj|2 008 ) 

The differential log N - log S relation are then integrated to- 
ward lower fluxes to produce the log N - log S relations in the 
integral form, A^(> S), where the number density of sources 
with fluxes above S is plotted (Figure 10). The achieved sen- 
sitivities, defined as the fluxes at which the sky area falls 
below 1% of the maximum area (1.14 deg 2 ), are 6 x 10~ 16 , 
-16 i v in-15 ,„h <; v m-i5 o,„ c-i 



be consistent with previous results with in uncertainties (e.g., 



8 x 10" lb , 3 x 10" 15 , and 5 x 10" 15 erg cm" 2 s" 1 in the 0.5- 
2 keV, 0.5-4.5 keV, 2-10 keV, and 4.5-10 keV band: above 
these fluxes (88 ± 5)%, (75 ± 3)%, (74 ± 4)%, and (52 ± 3)% 
of the XRB are resolved, when we adopt the XRB intensity of 
(7.5±0.4)x 10" 12 , (15.3 ±0.6) x 10" 12 , (20.2±l.l)x 10" 12 , 
and (12.3 ±0.7) x 1 0~ 12 erg cm" 2 s" 1 d eg" 2 , respectively, 
taken from Table 6 of ICarrera et al.l d2007l) with an appropri- 
ate flux conversion for the 4.5-10 keV band. For the esti- 
mate of the contribution of bright sources that is not well con- 
strained from the S XDS data, we have used the formula by 
ICarrera et al.l ((2007), who complied the results from Chandra 
deep surveys, XMM-Newton serendipitous surveys, and the 
ASCA medium sensitivity survey. 

Figure 1 1 shows the same log N - log S relations scaled 
by (5/10~ 14 ) 15 to stress the deviation from a Euclidean 
slope. The dashed curves represent the best-fit formula by 
ICarrera et al.l (12007); they are plotted in the figure covering 
only the region where the fit is performed. Our results are 
in good agreement with the previous work within «10%, 
indicating systematic errors, if any, is confined within this 
level . The red points are the results of the COSMOS sur- 
vey dCappelluti et al.l 12007b . Here we convert the result by 
I Carrera et al.l (120071) in the 4.5-7.5 keV band and that by 
Cappe lluti et al.1 (120071) in the 5-10 keV band to the 4.5-10 
ke V band assuming a photon index of 1 .5 (the flux conversion 
factor is 1.686 and 1.124, respectively). As seen from the fig- 
ure, the COSMOS source counts in the 2-10 k eV are 20-30% 
smaller than ours and the Car rera et al.l (120071) results. This is 
most probably be cause while our results a re truly based on the 
2-10 keV survey, Cappel luti et al.l (12007) make this plot from 
the 2^4.5 keV survey by converting the flux into the 2-10 
keV band assuming a photon index of 1.7. This could easily 
miss a population of sources with hard spectra or underesti- 
mate their fluxes. In the 0.5-2 keV band, our source counts at 
S < 3 x 10~ 14 erg cm " 2 s " 1 are systematically lar ger than both 
ICarrera et all (12007b and lCappelluti et all (120071) but are con- 
sistent with t he latest results from the XMM-Newton Lockman 
Hole survey dBrunner et al.ll2008l) . 

5. CLUSTERING OF SOURCES 
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(a) 0.5-2 keV 



0.03 

2.00 - 
1.00 - 



0.10 ; 
0.05 - 

0.03 L 




0.50 

CO 
QJ 

'-' 0.20 
to 

< o.io 

0.05 



10" 



io- 



Flux (ergs cm 2 s 1 ) 

1 1 1 I 1 1 " 1 1 ' ' I "" ' 
(c) 2-10 keV 



io- 15 io- 14 io- 

Flux (ergs cm- 2 s -1 ) 



(b) 0.5-4.5 keV 




10- 

Flux (ergs cm -2 

i 1 " 
(d) 4.5-10 keV 




io- 15 io- 14 io- 

Flux (ergs cm -3 s -1 ) 



FIG. 8.— The survey area given as a function of flux in the whole SXDS field in the (a) 0.5-2, (b) 0.5^.5, (c) 2-10, and (d) 4.5-10 keV band. 
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FIG. 9. — The differential log N - log S relations in the (a) 0.5-2, (b) 0.5^t.5, (c) 2-10, and (d) 4.5-10 keV band. The attached errors correspond to Icr Poisson 
error one in the number of sources. The solid lines are the best-fit broken power law model (see equation [T]and Table 7). 
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FIG. 10.— The log N - log S relation obtained from the whole SXDS field in the (a) 0.5-2, (b) 0.5—4.5, (c) 2-10, and (d) 4.5-10 keV band. 
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FIG. 11 .— The log N - log 5 relation obtained from the whole SXDS field in the (a) 0.5-2, (b) 0.5-4.5, (c) 2-10, and (d) 4.5-10 keV band, scaled by 
(5/10-' 4 ) 15 . The error bars in our data correspond to lcr statistical errors. The dashed curves correspond to the best-fit model by Carrera et al. 12007) obtained 
by using Chandra, XMM-Newton, and ASCA results. The open circles (red) are the results obtained from the COSMOS survey I Cappelluti et al. 2007). In (d), 
both results are converted into the 4.5-10 keV band assuming a photon index of 1.5. 
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FIG. 12. — The spatial distribution of detected X-ray sources in the whole 
SXDS field. Crosses, dots, and diagonal crosses correspond to those detected 
in the 0.5-2 keV, 0.5^1.5 keV, and 2-10 keV band, respectively. 

5.1. Spatial Distribution 

Investigating the evolution of the large scale structure of the 
universe is one of the main scientific objectives of the SXDS. 
From X-ray surveys, we can study the clustering properties of 
AGNs including both unobscured and obscured populations 
over a wide redshift range as a tracer of the large scale struc- 
ture. This also helps us to understand the AGN phenomena 
and their environment by measuring the mass of underlying 
dark matter halo where the AGN activity took place. 

Many studies have shown the presence of the signifi- 
cant large scale structure in the spatial (both 2-dimensional 
and 3-dimensional) distribution of X-ray sources (e.g., 
Vikhlinin and Forman 1995; Carrera et al. 1998; Akylas et al 
2000t lYang et all 120031; iGilli et al l 120031: iBasilakos et al 
20041 iMullis et al.l l2004: BasilakoseTalJ l2005b iGill 



2005b lYang et all l2006t iPucce tti 2006; iGandhi et alJ 120061: 



Carrera et alj 



et al 



2007; Miyaji et al. 2007). This leads pencil 



beam surveys to be inevitably subject to cosmic variance. The 
SXDS, thanks to its wide and continuous area coverage, pro- 
vides us with an ideal opportunity to investigate this issue. 
Figure 12 shows the location of sources detected in the 0.5-2, 
0.5^4.5, or 2-10 keV band (with ML>7) in our source cat- 
alog. Figure 13 shows the log N - log S relation separately 
derived from the seven pointings. As done in the analysis of 
§ 3.3, here we excluded any data with shorter exposures in 
overlapping regions of multiple pointings, and hence the re- 
sults are statistically independent one another. It can be seen 
that the source counts show significant variation among the 
pointings, depending on the flux limit. In particular, those of 
SDS7 in the 0.5-2 and 0.5^1.5 keV band are by more than 
30% smaller than the other fields in the medium to high flux 
range. These results indicate that the cosmic variance indeed 
exists over an area scale of one XMM-Newton FOV, ^0.2 
deg 2 , which is much larger than the FOV of Chandra obser- 
vatory. 

5.2. Auto Correlation Function 

The most straightforward approach to quantify the cluster- 
ing properties of sources is to calculate the angular two-point 
correlation function (or auto-correlation function; hereafter 
ACF). The ACF, w(9), is defined as 

dp = n 2 [i+ w (e)]dnidn 2 , 



where dP gives the probability of finding a pair of two ob- 
jects located at the differential solid angle of fli and O2 with 
a separation angle of 9, and n is the mean number density. 

To estimate the A CF, we essentially follow the same 
procedure as done in iMivaii et alj (12007b . We adopt the 
lLandv & Szalavl d 19931) estimator, 



w ohs (0) ■■ 



DD-2DR+RR 
RR ' 



where DD, DR, RR are the normalized number of data-data 
pairs, data-random pairs, and random-random pairs, respec- 
tively. In calculating the ACF, it is critical to ensure that the 
sensitivity map is correctly modeled at every position. Hence, 
considering possible position-dependent systematic errors in 
the sensitivity maps created in § 4.1, we introduce a "safety 
factor", /= 1.1, to multiply the nominal sensitivity limits 
and select only sources whose observed flux is higher than 
the (corrected) sensitivity limit at the detected position. The 
choice of / = 1.1 is quite conservative, and indeed we ob- 
tain essentially the same results even when / = 1 .0 is adopted. 
We make a random sample whose number is 10 times that of 
the actual data. The flux distribution of the random sample is 
taken to be the same as in the actual source list, and their posi- 
tions are randomly distributed by satisfying that the flux must 
be higher than the sensitivity limit at the allocated position. 

Figure 14 show our results of the ACF obtained in the 
0.5-2 and 2-10 keV bands. The flux range and number of 
sources used in the analysis are summarized in the 2nd and 
3rd columns of Table 8, respectively. The attached error is es- 
timated by simulation: we produce 100 sets of random sam- 
ples with the same source number as in the data, obtain a stan- 
dard deviation of the ACF signal obtained by the same proce- 
dure as applied for the actual data, and then multiplied it by 
VI + w(6). We fit the ACF with a power law form 

W (9) = (9/8 C ) 1 ^ 

in the range of 9 = 0.5 - 10 ', consi dering the size of the PSF 
of XMM-Newton. We fix 7 = 1.8 (lPeebleslll980l) . which is 
difficult to constrain by our data, and derive the correlation 
length 9 C as a free parameter. To take into account the cou- 
pling of the ACF between different bins, we utilize the "co- 
variance matrix", which is also obtained from the simulation 
of the random sample described above; we refer the read- 
ers to § 3.4 of IMivaii et al.l d2007l) for details. In the fit, we 
also consider an offset produced by the integral constraint, 
/ / w(9)d£l\dfl2 = 0, assuming that the power law form of 
the ACF holds over the whole area. The obtained values of 9 C 
are listed in the 6th column of Table 8. 

We detect an ACF signal from the 0.5-2 keV selected sam- 
ple with 9 C = 5.9"±o9». The signal is not significantly found 
from the 2-10 keV sample, however, with a 90% con fidence 
upper limit of 9 C < 1.5". By Limber's equation (Peebles 
1980), the ACF can be transformed to the 3-dimensional cor- 
relation function 



ar,z) = (r/r cfi r(l+z) 



-3-e+7 



once the redshift distribution of the sample is known. Here we 
assume the comoving clustering model where the correlation 
length r c ,0 is constant in the comoving coordinates, hence 
e = 7-3. Figure 15 shows the AGN redshift distribution for 
the 0.5-2 and 2- 10 keV samples estimated from the model by 
lUeda et all (120031) . From the 9 C values and redshift distribu- 
tion, we finally obtain r C) o = 14.9 ± 1.1 h~ l Mpc and r c .o < 7.6 
h~ l Mpc with a medium redshift of 1.3 and 1.1 for the 0.5-2 
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FIG. 13.— The log N - log 5 relation for each field in the (a) 0.5-2, (b) 0.5-4.5, (c) 2-10, and (d) 4.5-10 keV band. Black: SDS-1, Red: SDS-2, Blue: SDS-3, 
Green: SDS-4, Yellow: SDS-5, Cyan: SDS-6, Magenta: SDS-7. In the overlapping regions between different pointings data of lower exposure are excluded in 
the analysis i.e., the plots are statistically independent of one another. 



TABLE 8 

Clustering Properties of the X-ray Sources 



Sample Clustering Properties 



Detection Flux Range Number of $min-0max 7 e c "- d z cff b r c c - d 

Band (10~ 14 erg s" 1 cm 4 ) Sources (') (fixed) (") (/T 1 Mpc) 

0.5-2 keV 0.06-22 765 0.5-10 1.8 5.9^-2 1.3 14.9±1.1 

2-10 keV 0.32-33 573 0.5-10 1.8 0.1 « 1.5) 1.1 2.3 « 7.6) 

a The 2-dimensional correlation length, assuming a power law form of ( J^-)'~ 7 for the angular correlation func- 
tion. The fit is performed between 9 mm and # mii x with the integral constraint (see text). 
b The median redshift contributing to the angular correlation. 

c The correlation length assuming the comoving model (e = 7-3), converted from C via the Limber transforma- 
tion (see text). 

d The error is Icr while the upper limit is 90% confidence limit. 



and 2-10 keV selected sample, respectively. The bias param- 
eter, the ratio of the rms fluctuation amplitude rjg (Peebles 
1980) between AGNs and underlying mass, is found to be 
about 5.8±0.5 for the 0.5-2 keV selected AGNs at a median 
redshift of 1.3. 

The correlation length in the 0.5-2 keV is consistent with 
the previous results from s oft X-ray surveys at similar red- 
shift r ange r eported by e.g ., Basila kos et al.1 (120051) . iPuccettil 
(2006), and ICarrera et al.l d2007l) . although it i s larger than 
those obtained from the optically selected QSOs dCroom et al.l 
120051) . By contrast to several previous works that reported 
strong correlation of hard-band (>2 k eV) selected samples 
( Yan g et al.ll20"03t iBasilakos et aTl2004l) . we find no evidence 
for significant clustering signals in the hard band. The detec- 
tion of ACF in the 0.5-2 keV band but not in the 2-10 keV 
band is similar to what was found bv lGandhi et al.l d2006l) and 



ICarrera et al.l (120071) . It is curious that the first result from the 
XMM-COSMOS reports a relatively large correlation length 
for the 4.5-10 ke V band sample bu t smaller for the 2—4.5 
keV band sample dMivaii et all 120071) . Our 2-10 keV result 
is consistent with the COSMOS 2-4.5 keV result. We do not 
find significant ACF signals from the 4.5-10 keV sample, al- 
though the result is more subject to the statistical fluctuation 
due to the small number of sources detectable in this energy 
band. 

6. LONG TERM X-RAY VARIABILITY OF SOURCES 

We investigate the long term variability of X-ray sources 
detected in the SDS-4 field, for which two observations were 
performed with a time separation of about 2.5 year. The expo- 
sures of both observations are similar and are sufficiently long 
(> 20 ksec). Using the X-ray catalog produced from the com- 
bined images of the two observations, we examine the fluxes 
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FIG. 14. — The auto-correlation function of the SXDS sources in the 0.5-2 
keV and 2-10 keV bands. The errors are 1 <r. The line represents the best-fit 
model (power law plus constant to account for the integration constraint). 
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FIG. 16. — The comparison of the 0.5^1.5 keV flux measured in the two 
different epochs for sources detected in the SDS-4 field. 



0.8 



0.6 - 



N 



x) 0.4 - ' 




o.o r ^ 

1 2 3 4 5 

z 

FIG. 15. — The normalized redshift distribution of AGN in our samples 
selected in the 0.5-2 keV (s olid curve) and 2- 10 keV (dashed curve) bands, 
estimated from the model by Ueda et al. I 2 0"03h . 

of the cataloged sources separately in each epoch by fixing 
their positions. To ensure high signal-to-noise ratio, we only 
use sources whose summed ML value from the 0.5-2 keV 
and 2-4.5 keV bands exceeds 15 in the combined catalog and 
vignetting-corrected exposure is longer than 15 ksec in both 
observations. Figure 16 shows the comparison of the 0.5-4.5 
keV flux between the first and second epochs for this sam- 
ple. As noticed from the figure, many sources at intermediate 
fluxes that have sufficiently small error bars show a signifi- 
cant variability, demonstrating the impact of time variability 
in studying AGN properties. The results are consistent w ith 
previous studies dPaolillo et al.ll2004t iMateos et alll2007l) re- 
porting that the fraction of variable sources on time scale of 
months to years is 80-90% or higher in Chandra and XMM- 
Newton deep fields. 



7. SUMMARY 

We have presented the first X-ray source catalog and their 
basic X-ray properties in the SXDS field, based on the seven 
XMM-Newton pointings performed between 2000 and 2003 
that cover a continuous area of 1.14 deg 2 in the 0.2-10 keV 
band. The catalog consists of 1245 sources in total, consist- 
ing of those detected either in the 0.3-0.5 keV, 0.5-2 keV, 
2^.5 keV, 4.5-10 keV, 0.5-4.5 keV, and 2-10 keV with de- 
tection likelihood larger than 7. From empirical sensitivity 
maps derived by simplified simulation, we obtain log N - log 
S relations in various energy bands bands, which bridge the 
flux range between Chandra deep surveys and brighter sur- 
veys. Clustering properties of X-ray sources are also studied 
by means of auto correlation function. We detect significant 
signals that can be fit by (6/0 c T S with C = 5.9'% f, for the 
0.5-2 keV selected sample at an estimated median redshift of 
Zett « 1.3, and an 90% upper limit of 6 C < 1.5" for the 2-10 
keV selected sample at z e g « 1.1. Two pointing data sepa- 
rated by 2.5 years indicate the importance of flux variability 
on a time scale of years. Our results establish the average 
properties of X-ray sources at fluxes where a majority of the 
XRB is produced, being least affected by cosmic variance. 
The data presented here constitute a major component of the 
SXDS project and shall be used for a number of research pro- 
grams in combination with other wavelengths data. 
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TABLE 2 

The X-ray source list of the SXDS 



Position (J2000) Likelihood Count Rate Hardness Ratio 



# 


R.A. 


Dec. 


7/ 


0.3-0.5 


0.5-2 


2-4.5 


4.5-10 


0.5-4.5 


2-10 


0.3-0.5 


0.5-2 


2-4.5 


4.5-10 




HR2 


HR3 


Field 


Offset 


Exp. 


Bgcl. 


Note 


























-1) 




















0001 


33.85387 


-4.92870 


2.36 


0.0 


19.9 


3.2 


0.4 


26.5 


5.1 


0.00± 0.09 


3.45± 0.79 


1.39± 0.63 


0.86± 0.81 


1.00±0.05 


-0.42±0.21 


-0.24±0.49 


5 


15.6 


9.6 


2.6 




0002 


33.85392 


-4.90221 


2.24 


16.1 


33.7 


2.5 


0.8 


37.1 


6.2 


2.69± 0.74 


6.33 ± 1.25 


1.38± 0.72 


1.18± 0.90 


0.40±0.14 


-0.64±0.16 


-0.08±0.46 


5 


16.0 


5.5 


2.6 




0003 


33.86346 


-5.03354 


2.88 


0.3 


0.0 


0.0 


6.8 


0.4 


7.5 


0.24± 0.27 


0.18± 0.31 


0.09 ± 0.25 


2.72± 0.82 


-0.13±1.01 


-0.33±1.43 


0.94±0.17 


5 


14.8 


12.2 


2.6 




0004 


33.87124 


-4.96292 


2.99 


5.2 


14.1 


0.0 


0.0 


11.1 


0.0 


0.79± 0.34 


2.33± 0.57 


0.00± 0.11 


0.00± 0.28 


0.49±0.19 


-1.00±0.09 




5 


14.2 


15.5 


2.6 




0005 


33.87257 


-4.89099 


1.85 


2.5 


11.8 


6.7 


0.7 


22.4 


8.6 


0.68 ± 0.38 


2.85± 0.73 


1.76± 0.55 


0.88 ± 0.69 


0.61±0.19 


-0.24±0.19 


-0.33 ±0.38 


5 


15.2 


9.7 


2.6 




0006 


33.87878 


-4.98760 


2.84 


0.0 


2.6 


2.9 


0.0 


8.7 


2.8 


0.00± 0.04 


0.89± 0.41 


0.61± 0.29 


0.00± 0.38 


1.00±0.10 


-0.19±0.32 


-1.00±1.24 


5 


13.6 


17.9 


2.5 




0007 


33.89418 


-5.06675 


0.56 


131.3 


554.7 


5.5 


0.0 


484.2 


4.3 


6.86 ± 0.88 


25.08 ± 1.56 


1.20± 0.43 


0.03 ± 0.25 


0.57±0.05 


-0.91 ±0.03 


-0.94±0.40 


5 


13.6 


14.1 


2.8 




0008 


33.90476 


-4.94904 


1.25 


24.9 


41.6 


6.7 


0.0 


50.9 


5.9 


1.56± 0.33 


3.57 ± 0.56 


0.85 ± 0.27 


0.00± 0.21 


0.39±0.11 


-0.62±0.11 


-0.99±0.49 


5 


12.3 


16.8 


2.2 




0009 


33.90754 


-5.08932 


2.64 


0.0 


6.8 


3.4 


0.0 


13.6 


3.0 


0.00± 0.06 


1.75± 0.56 


0.77± 0.34 


0.14±0.32 


1.00±0.06 


-0.39±0.23 


-0.69±0.61 


5 


13.4 


17.3 


2.7 




0010 


33.91171 


-5.01450 


1.30 


10.6 


72.1 


6.0 


2.3 


75.5 


11.0 


0.89± 0.25 


4.53 ± 0.59 


0.73± 0.25 


0.90± 0.44 


0.67±0.09 


-0.72±0.09 


0.11±0.29 


5 


11.8 


17.4 


2.2 





The pointing ID (1-7) where the source is detected. 

The offset angle in units of arcminutes from the mean optical axis in the corresponding pointing. 

The total pn-equivalent exposure (sum of pn, MOS1, and MOS2) at the source position in units of ksec, corrected for the vignetting in the 0.5^1.5 keV band. 
The background rate in the 0.5^1.5 keV band at the source position in units of 1CT 3 counts ksec -1 arcsec~ 2 . 

Set to be "C" if multiple-source fitting was performed to obtain the fluxes to take into account possible source confusion with nearby sources (see § 3.2.3). 



